Probabilistic Phase Space Trajectory Description for Anomalous Polymer Dynamics 
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It has been recently shown that the phase space trajectories for the anomalous dynamics of a 
tagged monomer of a polymer — for single polymeric systems such as phantom Rouse, self-avoiding 
Rouse, Zimm, reptation, and translocation through a narrow pore in a membrane; as well as for 
many-polymeric system such as polymer melts in the entangled regime — is robustly described by 
the Generalized Langevin Equation (GLE). Here I show that the probability distribution of phase 
space trajectories for all these classical anomalous dynamics for single polymers is that of a fractional 
Brownian motion (fBm), while the dynamics for polymer melts between the entangled regime and 
the eventual diffusive regime exhibits small, but systematic deviations from that of a fBm. 
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In its terminal relaxation time t, t ^ N'^^ a polymer of 
length N displaces itself in space by its own size, which it- 
self scales as ~ N^ [1, 2]. The values of k and ^ vary from 
system to system. E.g., for phantom (self-intersecting) 
polymers ^ = 1/2, and for self-avoiding polymers ^ = v, 
with v = 3/4 in two and w 0.588 in three dimensions 
respectively. Similarly, for polymer dynamics in the ab- 
sence of hydrodynamic interactions (Rouse: k = 1 -(- 2^ 
3, |j|), and polymer dynamics in a good solvent (Zimm 
l|,l3,l5|: K — 3^). The above means, in the simplest case, 
that the mean-square displacement (MSD) of a tagged 
monomer of a polymer must behave ^ ^^^/k ^j^f^[\ time 
T, and ^ t thereafter. Since 2^/k is not necessarily 
unity, the dynamics of a tagged monomer in a polymer 
is anomalous till time r, and the polymer's diffusion co- 
efficient scales ^ N'^^~'^. 

Starting from a microscopic description, there are 
two main approaches to model anomalous dynamics in 
stochastic systems [6]: (i) Continuous Time Random 
Walk (CTRW) and the associated fractional Fokker- 
Planck Equation (fFFE), providing a probabilistic de- 
scription of phase-space trajectories, and (ii) the Gen- 
eralized Langevin Equation (GLE) [8|, which describes 
individual phase space trajectories. Physical systems ex- 
hibiting anomalous dynamics, for which probabilistic de- 
scription of phase-space trajectories as well as description 
of individual trajectories can be obtained, are not only 
relatively rare, but also, relating one description to the 
other often requires approximations [9|. Anyhow, given 
the ubiquity of anomalous dynamics in polymeric systems 
[l|, y|, one would expect them to have been thoroughly 
examined from this perspective. To the best of my knowl- 
edge however, probabilistic description of phase space 
trajectories for anomalous polymer dynamics has only 
been considered for isolated cases such as Rouse chain 
[lCr|, and polymer translocation |lll4l5| : in [ill Il2|. an 
fFPE approach has been put forward, wherein, a (power- 
law) waiting time before each move of the monomer in 
the pore is assumed to cause the anomalous dynamics. 



This approach is at odds with numerical studies by oth- 
ers [13|, who report that, for the translocation of an in- 
finite polymer, the probability distribution is Gaussian 
in space, but with a width that scales anomalously in 
time. Further, tire anomalous dynamics of translocation 
has been shown to match that of the fractional Brownian 
motion (fBm) lj,ll5|, which is also in contradiction with 



the fFPE approach. 

As for the description of individual trajectories for 
anomalous polymer dynamics, two recent papers |l6l.ll7| 
show that without external forces, the motion of a tagged 
monomer of a polymer — for phantom and self-avoiding 
Rouse, Zimm, reptation, translocation, and polymer 
melts — is robustly described by a unified Generalized 
Langevin Equation (GLE). The force (j){t) experienced by 
the tagged monomer is related to its velocity v{t) via 



m 



dt'fi{t-t')v{t')+g{t), 



(1) 



where /i(t) is the memory kernel, and the noise g{t) sat- 
isfies {g{t))a = and the fluctuation-dissipation the- 
orem (FDT) (g.(i)5A(i'))o = ksTd^xfiit ~ t'), with 
a, X — (x, y, z). Here (. . .)o denotes an average over the 
noise realizations, including an average over equilibrium 
configurations of the polymers at i = 0~. Further, with 
v{t) responding to (p{t) as 'yv{t) ~ (pit) + /(i), where 7 
is the (effective) damping coefficient for monomeric mo- 
tion, and f{t) is a random force satisfying (/(t)) = and 
the FDT {Uit)h(t')) = 2jkBTS^xS{t - t'), one has 



vit) = 7-1 



dt'^it - t')v{t') + g{t) + fit) 



(2) 



In this formulation, /^(i)~ t "e *''^ for some < a < 1 
(for a list, see Table I of 0). The FDT then ensures that 
the MSD ~ t" till time r, and -- t thereafter. This for- 
mulation also robustly yields the correct drift behavior of 
a tagged monomer under weak external forces [17|, like 
the Nernst- Einstein relation. The GLE (HE]), describ- 



polymeric system 


derivation of Eq. (g]) [2] 
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phantom Rouse 


directly from Rouse equation 
no approximation 


2 


70 = N'y 
7p^o = 2Af7 


Gix'^kBTp^/N 


Zimm (phantom) and 
polymer in a ^-solvent 


Smoluchowski equation; pre-averaging 
approximation on the mobility matrix 


3/2 


70 = 3(6^^iV)i/2j,^/8 
7p^o = {Un^Np^^^s 


6tt^ keTp'^/N 


Zimm (self-avoiding) 


Smoluchowski equation; pre-averaging 
approximation on the mobility matrix 


3// 


70 = rjsN" 
7p = VsN-'p'-' 


N-'^''kBTp^+'"' 


reptation 
(curvilinear co-ordinate) 


entropic (curvilinear) chain tension; 
i.e., Eq. H is only one-dimensional 


2 


70 = iV7 
7p^o = 2N-f 


Gix^kBTp^/N 



TABLE I: Elaboration of Eq. ([4]) for phantom Rouse, Zimm, polymer in a ^-solvent and reptation. Here r]s is the solvent 
viscosity. Where applicable, the effective damping coefficient 7 for monomeric motion, as defined in Eq. ([2]), is oc 773. Note, in 
all cases, that k, is simply the power of p in the expression kp/-^p for p 7^ 0. 



ing non-Markovian trajectories in phase space, demon- 
strate that there is no power-law waiting time (assumed 
in the modehng of translocation by the fFPE); instead, 
the anomalous dynamics stems from the fact that each 
move of the tagged monomer tends to be undone later. 

Here I demonstrate that, with a (5- function distribution 
in 3D space at i = 0, the probability distribution of a 
tagged monomer for phantom Rouse, self-avoiding Rouse, 
Zimm, polymers in a 0-solvent, reptation, translocation, 
and polymer melts in the entangled regime is given by 



P(r,t|ro,0) = e-(*'-*''')'/2A(t)/j27rA(t)]3/2, 



(3) 



with A(t) = At°e~*/^ -f BN'^i-'^t, where A and B are 
two system parameters-dependent constants, and N'^^~'^ 
is the scaling of the polymer's diffusion coefficient with N 
[for translocation, Eq. ^ holds only until the polymer 
disengages from the pore; see later]; i.e., the anomalous 
dynamics of an infinite polymer is that of the fBni [18]. 
One should keep in mind that Eq. ([3|) is demonstrated 
here for polymeric systems wherein the polymers are far 
away from any boundary (including translocation), as 
Eq. ^ cannot be reconciled with nontrivial boundary 
conditions. The CLE ^HQ^ and Eq. ^ thus provide 
a complete description (i.e., of individual trajectories as 
well as that of trajectory distribution in phase space) of 
anomalous dynamics in polymeric systems; as mentioned 
earlier, this is relatively rare for physical systems. 

In fact, P(r, t|ro,0) for phantom Rouse, Zimm, poly- 
mers in a ^-solvent, and reptation can be obtained analyt- 
ically, irrespective of the GLE description, thanks to the 
fact that their dynamics is described by that of the poly- 
mer's fluctuation modes. For these systems, the location 
Rn(t) of monomer n can be expressed, in terms of the 
mode amplitude {Xp{t)}s for p = 0, 1, 2 . . ., as Rn{t) = 
^o{t) + S^i -^p(^) cos{npn/N), obeying the boundary 
condition that the polymer's chain tension vanishes at the 
free ends, i.e., [dRn{t)/dTi]n=o — [dRn{t)/dn]n=N — 0. 
As the {Xp{t)}s can be obtained by an inverse cosine 
transformation of the above, the polymer dynamics is 



simply reconstructed from the LE satisfied by each spa- 
tial component (denoted by a) of the {Xp{t)}s ^ 



lpXpa{t) = -kpXp^it) + fpait), 



(4) 



where the stochastic force satisfies (fpait)) = and the 
FDT ifpafqx) = 2jpkBT5ax6pgS{t-t'). The list of 7^ and 
Kp values for these systems appear in Table HI note that 
the relaxation time Tp for the p-th mode Tp = jp/kp — 
(,-ip-K £qj. j, ^ 0; one can calculate c from the table. 



polymeric system 


A(f) 


phantom Rouse 


~ t^/^ till r and ~ t thereafter 


Zimm (phantom) and 
polymer in a ^-solvent 


~ t^/^ till r and ~ t thereafter 


Zimm (self-avoiding) 


~ t^/^ till r and ~ t thereafter 


reptation 
(curvilinear co-ordinate) 


~ t^'^ till r and ~ t thereafter 


self-avoiding Rouse 


~ 1 1+2'^ till r and ~ t thereafter 



TABLE IL Systems of TableUand their A(t)-behavior. Note 
that the scaling exponent of A(f) with time, in each case, is 
given by 2^/k, as noted in the first paragraph of the paper. 

With the above, P{r,t\ro,0) is calculated as follows. 
For fixed {Xp }s that correspond to the position tq of 
the tagged monomer at t = 0, one determines the prob- 
ability 7'(Xp,t|X^°\ 0) from Eq. (HI). One then obtains 
P{f, t\ro, 0) by integrating over all values of {Xp{t)} that 
correspond to the position r of the tagged monomer at 
time i, taking into account the equilibrium distribution of 
the {Xp }-values. With the tagged monomer being the 
middle monomer, this calculation, demonstrating Eq. ([3]) 
for phantom Rouse, Zimm, polymers in a 0-solvent, and 
reptation are detailed in Appendix A: therein it is seen 
that [c.f. Eqs. (|A13IIA17p ] 



A{t) 



dp 



-t/rf 



(5) 



i.e., the exponent of A(i) in time is determined by the 
integral in Eq. ([5]). The corresponding scalings of A{t) 
summarized in Table HIl Further, based on [19], I note 
that one can construct an effective Eq. (U) also for self- 
avoiding polymers, with p-independent 70 and "fp^o, and 



altering the polymer's chain tension locally [Eq. ^]. 
In response to this strain, in subsequent times there 
is an enhanced chance, for the monomer, to undo the 
move [Eq. ^]. This physics is best represented by 
discretizing the movement of the middle monomer in 



ti) for some m, and io = by choice. In between 
these moves the middle monomer remains stationary, 
i.e.. the dynamics of the polymer is given by those 
of the mode amplitudes {Yp{t)}s, with the monomer 
locations r„(i) relative to the middle monomer ex- 



pressed as fn{t) = 4 ^ Fj 



i?W^ 



4Ep?p''(i)sm2M)lie(7V/2 



kp ~p^+^''. The corresponding scaling behavior for A{t), time for the ensemble £ [17|; e.g., v{t) — ^^QSfiS{t 
also calculated in Appendix A, is listed in Table [TTl and 
is verified by simulations in Fig. [T] 

It is interesting to note here that the mode expansion 
technique to establish Eq. ([3]) for the systems of Table 
|T] implies that Eq. ([31) is a simple consequence of the 
physical connectivity of the polymer chain. 

The expansion of the monomer co-ordinates into poly- 
mer's fluctuation modes also provides an insight into how, 
given the GLE (ll][2|), one may expect P(r, i|ro,0) to be 
Gaussian. E.g., consider an ensemble £ of polymers, at 
equilibrium at t = 0~, with a given velocity history of 
the middle monomer between times and t. For this en- 
semble, if [g{t) + f{t)] can be shown to be Gaussian, then 
the displacements of the middle monomer in an infinites- 
imal time between t and {t + dt) is Gaussian distributed 
about the mean — 7^^ /„ dt'ij,{t—t')v{t'), since v{t) is pro- 
portional to [g(t) + f(t)]. Such (infinitesimal) Gaussian 
displacements, accumulated over time, would then mean 
that P{rcr,t\rcro,Q) has to be Gaussian. 
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'(i)sin^i^^iii^e(n-A^/2)- 

-j^ v.^.,/^ n) for p = 0,1,2,... 

(the superscripts for Y{t) correspond to the right and 
the left halves of the polymer). The {Y^(t)}s are ob- 
tained from {rn(i)}s via the inverse sine transform, and 
are readily shown to satisfy the boundary condition that 
the chain tension vanishes at the open ends of the poly- 
mer. Through this formulation, the polymer's chain ten- 
sion at the middle monomer, expressed in terms of the 
{Yp{ty\s, changes discretely at {ii}s, while in between, 
the relaxation of the chain tension gives rise to the mem- 
ory kernel of Eq. ([T|) . To work this out for all systems of 
Table HI one needs to re-perform, as applicable, the pre- 
averaging approximation in terms of the {Yp{t)}s, which 
is a cumbersome task. For the sake of simplicity, I there- 
fore only consider the phantom Rouse case here; for this 
system, in times U < t < tj+i, the {Yp(t)}s for each half 



independently obey the LE 17 1 
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(6) 
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FIG. 1: (color online) Equation ([3]) for self-avoiding Rouse 
polymers (A^ = 400): model details can be found in Appendix 
B. Left: second moment m2{t) = '^^{[ra{t) — rCT(0)]^)/3 

of the distribution ((3}, scales as 1^+^" ~ t'^'^'^ till r, and 
~ t thereafter. Right: m2(t) against raiit) = '^^{[ra{t) — 
ra(0)]*)/3. Red line: m2(t) = [m4(t)/3]^/^ — note that 
mA{t) = 2>[ni2{t)\ for a Gaussian distribution. Data obtained 
from a time-series of 200,000 consecutive snapshots, separated 
by 400 time units each, of 256 different polymers. 

In order to mathematically appreciate the Gaussian 
behavior of [g{t)+f{t)\, one needs to recall the physics be- 
hind the anomalous dynamics for these systems [IGJ, [l7 1 : 
a move of the middle monomer creates a local strain by 



w_ith 7p = 27V7, qp = 6TT^kBT{2p+lf/N, T2p+i = -fp/qp, 
(hp) = and {hpa{t)hqx{t')) = ^pkBT5{t - t')5pg5ax. 
Then g{t) in Eq. dTP is given by (see Eq. (19) of [17]) 

5W =4^^^(?^|e-/--M>;<'-H0-) + r«(0-)] 
p I 

+ i f dt' e~'^'-''^'^'-^^[h^;\t') + h^^\t')]\. (7) 



IpJa 

With /(t), {/ip(t)}s, {yj,(0~)}s being Gaussian dis- 
tributed with zero mean [as of Eqs. © and ^ respec- 
tively], [g{t) -f f{t)] also has to be Gaussian. 

Thus, to summarize so far, having expanded the 
monomer co-ordinates in polymer's fluctuation modes, I 
have shown that P(r, i|ro,0) for phantom Rouse, Zimm, 
polymers in a 0-solvent, reptation and self-avoiding 
Rouse is Gaussian; and have illustrated, for the specific 
case of phantom Rouse, that Gaussianity of P(r, i|ro,0) 
is expected from the GLE ([l][2)) since the noise term 
[g{t) + f{t)] is Gaussian. Unfortunately however, the 
mode expansion does not work as simply for transloca- 
tion (for which, the tagged monomer is the one within 



4 10^ 10^ 10^ 10 

10 hi iiiiiiii I iiiiiiii I iiiiiiii I iiiiiiii I iii^i 



(a) 



10^ 



10 



10 




io' lo' lo' lo' lo' lo' «4(^) 
t 



10 



13 5 7 9 

10 10 10 10 10 

HIIIIII IIIIIIII IIIIIIII IIIIIIII IIIIIIII IIIIIIII IIIIIIII lit 




I I I I I I I III 



lo' lo' lO' 10^ '«4(^) 
t 



FIG. 2: (color online) (a) Left: The second moment m2{t) — (1/3) ^^{[ra-{t) — rcr{0)]^) of the distribution ((3| for the middle 
monomer of a tagged polymer in a melt. Right: 7712 (i) against m4(i) = {l/3)^^{[ra{t)—rtT{0)]''^)',iedlme: m2{t) — [m4(i)/3]^'^. 
Data obtained from long time-series of 155,000 consecutive snapshots, separated by 10* time units each, of 1,728 different 
polymers with A*' = 1500. See text for graph description, (b) Left: The second moment ni2{t) — [s{t) ~ N/2]^, where s{t) 
is the monomer number threaded in the pore at time t, s(0) — N/2. Right: m2{t) against m4{t) — {[s{t) — N/2]'^); red line: 
m2{t) = [m4{t)/3]^'^ . Data averaged over 8,192 different polymers of A'^ = 1000, for the turquoise points in the left graph at 
least one polymer of the 8,192 polymers has disengaged from the pore. See text for graph descriptions. 



the pore at time t, i.e., it does not even have a fixed in- 
dex) and for polymer melt. Neither can Gaussianity of 
[g{t) + f{t)] be shown, and therefore, one has to rely on 
computer simulations. The data for the melt are pre- 
sented in Fig. HJa): simulations are performed at overall 
monomer density unity; the simulation details, same as 
that of [17| where the GLE has been shown to describe 
the dynamics of the middle nronomer of a tagged poly- 
mer in the entangled regime, can be found in Appendix 
B. Following reptation theory for polymer melts, one ex- 
pects the second moment m2(t) to behave ~ t^/^ in the 
entangled regime, which starts around time ~ 10^ for this 
model; in this regime an effective exponent 0.27 is found 
(the red line in the left graph). While the right graph is 
consistent with Eq. ([3]) in the entangled regime, betweeir 
the entangled regime and the eventual diffusive regime, 
P{f, t\fo, 0) does deviate very slightly, and systematically, 
from ([3]), indicating that the noise term [g{t) + f{t)] is 
not Gaussian during this time. 

The data for unbiased translocation are presented 
in Fig. m^b) [see model details in Appendix B]. The 
monomer number within the pore at time t is denoted 
by s{t). Polymers are equilibrated with s(0) — N/2. 
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Here m2{t) ^ ^(i-i-i')(i+2i/) ^ ^0-73^ ^^^ should cross over 
to diffusive behavior, as predicted in [20j, although the 
crossover is slow: at long times polymers start disen- 
gaging from the pore, corresponding to the flattening 
of m2{t) (the turquoise points in the left graph), hence 
the true diffusive behavior can only be observed for very 
long polymers. Nevertheless, Eq. ([3]) is verified cleanly 
up to the point when all polymers remain threaded in 
the pore: this behavior is in agreement with [13, 14] [al- 
beit they report an anomalous exponeirt different fronr 
(1 + i^)/(l + 2i/)]. and contradicts [13[, reaffirming that 
fFFE is not applicable for polymer translocation. 

Finally, 1 irote that the drift of a tagged monomer due 
to a weak external force F (i.e., in the linear response 
regime), is given by f^F till time r and tF thereafter, 
where a is the anonralous exponent. Describing this re- 
quires a simple extension of Eq. ([3]) . 

Ample computer time from the Dutch national super- 
computer facility SARA, and help from Gerard Barkema 
with simulations are gratefully acknowledged. It is also 
a pleasure to thank Gerard Barkema and Rainer Klages 
for helpful comments on the manuscript. 



Appendix A: : Derivation of Eq. (I3|) for the middle monomer for phantom Rouse, Zimm, polymers in a 

^-solvent, reptation, and self-avoiding Rouse polymers 



I tag the nriddle nronomer of the polymer, and here I obtain Eq. ^ for its dynamics. 



I start with the Langevin equation ^ describing the evolution of the p-th mode ampUtude (p = 0, 1, 2, . . .), viz., 

(Al) 



Ip 



— kpJypa-[t) + Jpcr{t) 



dt 



with (/p^) = and the FDT {fpa{t)fq\{t')) = 2^pkBTS{t-t')6pqSa\. As noted in TableU Eq. ^M^ can be derived for 
phantom Rouse, Zimm, polymers in a ^-solvent and reptation. A straightforward result that follows from Eq. (jAip 
is that 



{Xp4t)Xgx{t')) = 6pg6ax{kBT/kp)e-''''^'-'">^''- = SpqS^),{kBT/kp)e-^'-''^^^- 



(A2) 



where Tp is the relaxation time of the p-th mode {p ^ 0) for the polymer. When Eq. (jA2[l is combined with the 
corresponding time correlation function for mode amplitudes for self-avoiding polymers ^19l] . namely {Xp{t) ■ Xq{t')) oc 
]\j-2up-{i+2iy)^-t/Tpg^^ with Tp ^ {N/pY^"^^ , one can formulate an effective Eq. (jj]) with both 70 and ^p^Q independent 
of p, and kp ~ p-i'^+'^'^) [this implies that for a self avoiding polymer Tp ^ (N/p)^^'^'^]. Given this, I will henceforth 
use Eq. (jAip also for self-avoiding Rouse polymers. 



The Fokker- Planck equation for the probability V(Xpa,t) that corresponds to the LE is given by [22 1 



dV{Xp,,t) 
dt 



Ip 



d 



dXpa- 



[Xp,r{Xp,,t)] + 



kgT d^r{Xp,,t) 



=Tp ^ for p#o 



7p 
=aj,Tp^ for p^O 



dXl, 



(A3) 



where Op = ksT/kp for p ^ 0, and oq = ksT/jQ. The solution of Eq. (|A3I) . with the initial condition that 
P{Xpcr, 0) = S{Xpcr — XpJ\ is obtained as follows. 



(i) For p = 0, fcp = i.e., (|A3p is a simple diffusion equation. Its solution is given by [21] 

1 



PiXoa,t) = 



y/27raot 



exp 



{Xpa - Xq^ ) 

2aot 



(A4) 



(ii) For p 7^ 0, Eq. 



2j) can be verified by direct substitution of its solution 

{Xp,~X^"Je-'/-^r 



V{Xp„,t)^ 



1 

27rap(l-e-2*/^p) 



exp 



2ap(l-e-2*/rp) 



(A5) 



Next, as noted above Eq. (|1]), in terms of the mode amplitudes, the location of the middle monomer (n = A^/2) at 
any time t is given by 



r{t) ^ Xp{t) + 2Y,Xp{t) cos ^ 
p=i 



Using Eq. (|A6p , I obtain, upon averaging over all possible initial states of the polymer at i = 



P(r.,i|ro., 0) = n / dX^'i^ 'P.^i^'^J) S 

p=0"'-°° 



roa 



^^°^ + 2E^i°^-°« 



q=l 



(0)rr..^ 



/oo " / 00 \ ' 

dXp„ ViXp^,t) 6 r^- iXoa + 2j2 Xqa cos y , 



(A6) 



(A7) 



where Pcq(A") is the equilibrium probability of X, i.e., a Gaussian, obtained by taking the i — > 00 limit of Eq. (jA5 
At this stage, because of the (5-functions in Eq. (|A7[) . it is easiest to Fourier transform P(ro.,i|r-0(T, 0), defined as 



Vu,k',t = 2^ydr,dro.e^['-''^'+'=^-]pK,t|ro.,0], 



(A8) 



which reduces Eq. (jA7[) to 

r\^ —OO 

/oo 
-oo 



X 



(A9) 



At this point, in order to follow through the calculation of Vk.k'-.t, I need the two following integrals: 
(a) 



dXocr 



1 



VSTTOot 



exp 



{Xqo- - Xq^ ) 

2aot 



ik'X, 



Oct 



g-ixZ^k'-^aotk'^_ 



(AlO) 



(b) for p^O: 



dXr. 



cxp 



2ap(l-e-2*/-.) ^^'^^P-cos 2 



pvr 



27rap(l-e-2*/Tp) 

_ -2iX(^'e-*/^Pfc'cos(pir/2)-2ap(l-e"^*/^P)fc'2cos^(pir/2) 



(All) 



Using (a-b), I now integrate over X^J (i.e., the location the center-of-mass of the polymer) with a uniform probability 
density measure yields \/2TTS{k + k'), which leads me to 

oo r „oo 



p=i 






(A12) 



'Pk,k';t oc (5(fc + A;') implies that P{r^,t\roa,0) is a function of {r„ — roa)- 

Finally, I now need to evaluate the discrete sum in the exponent of Eq. (|Al2p . Having noticed that cos(g7r/2) = 
for odd g- values and cos2(g7r/2) = 1 for even (/-values, the sum can be converted into an integral; thereafter the 
inverse Fourier transform from k to {tc — roo-) leads to Eq. Q, with the behavior of A(i) presented in Table HIl 
With the corresponding scaling of jp and Tp = "fp/kp for phantom Rouse, Zimm, polymers in a 0-solvent, reptation, 
and self-avoiding Rouse polymers (see Table |l|, these integrals are listed below. Note that in Eqs. (|A13IIA17|) I omit 
constants in converting the discrete sums to integrals. 



A. Phantom Rouse: 



oo />oo ; 

Ay aqC0S^{qn/2){l - e-'^^^) ^ kgT / ^ {1 

B. Phantom Zimm and polymers in a ^-solvent: 

4 V a, cos2(qV2)(l - e-'/--) -^ ksT 4(1- 

9=1 -^0 T 

C. (self- avoiding) Zimm: 

oo /"OO 7 

4^a,cos2(W2)(l - e-*/-^) ^ ksT -^ (1 



-cq^t 



)^Vt. 



(A13) 



3-=?''^'*) ^ i2/3_ 



(A14) 



_ e-^^i *\ ^ y:2/3 



')^t' 



(A15) 



D. reptation (curvilinear co-ordinate): 

E. self-avoiding Rouse: 
4^ a,cos2(qV2)(l - e-*/"') ^ feBTJ ^ (1 - e 



;-^«'*) ^ Vt. 



9=1 



(A16) 



-cqi+^''t\ ^ ^2i//(l+2i/) 



(A17) 



Clearly, these power-law behavior of A(t) cannot hold longer than time r, this is also noted in Table HIl 

I 



Appendix B: Simulation details 

Over the past years, a highly efficient simulation ap- 
proach to polymer dynamics has been developed in our 
group. This is made possible via a lattice polymer model, 
based on Rubinstein's repton model ^] for a single 
reptating polymer, with the addition of sideways moves 
(Rouse dynamics). A detailed description of this model, 
its computationally efficient implementation and a study 
of some of its properties and applications can be found 
in|24|]. 

In this model, each polymer is represented by a sequen- 
tial string of monomers, living on a face-centered-cubic 
lattice with periodic boundary conditions in all three 
spatial directions. Hydrodynamic interactions between 
the monomers are not taken into account in this model. 
Monomers adjacent in the string are located either in 
the same, or in neighboring lattice sites. The polymers 
are self-avoiding: multiple occupation of lattice sites is 
not allowed, except for a set of adjacent monomers. The 
number of stored lengths within any given lattice site is 
one less than the number of monomers occupying that 
site. The polymers move through a sequence of random 
single-monomer hops to neighboring lattice sites. These 
hops can be along the contour of the polymer, thus explic- 
itly providing reptation dynamics. They can also change 
the contour "sideways", providing Rouse dynamics. Each 
kind of movement is attempted with a statistical rate of 
unity, which defines the unit of time. This model has 
been used before to simulate the diffusion and exchange 
of polymers in an equilibrated layer of adsorbed poly- 
mers [25|, dynamics self-avoiding Rouse polymers [26], 
polymer translocation under a variety of circumstances 
[la, I20I, l27| , and the dynamics of polymer adsorption [28| . 

The same model has been used for the polymer 
melt simulations (here the polymers are both self- and 
mutually-avoiding) for a system of size 60'^ with an over- 
all monomer density unity per lattice site. Due to the 
possibility that adjacent monomers belonging to the same 
polymer can occupy the same site, overall approximately 
40% of the sites typically remain empty. 

Initial thermalizations were performed as follows: com- 



pletely crumpled up polymers are placed in lattice sites 
at random. The system is then brought to equilibrium 
by letting it evolve up to 10^ units of time, with a com- 
bination of random intermediate redistribution of stored 
lengths within each polymer. Additional details on the 
melt simulations can be found in [17i |. 



[3] 
[4] 
[5] 
[6] 

[7] 



[1] De Gennes P-G, 1985 Scaling concepts in polymer physics 

(Ithaca, Cornell, revised edition). 
[2] Doi M and Edwards S F, 2003 The theory of polymer 

dynamics (Clarendon, Oxford) 

Rouse P E, 1953 J. Chem. Phys. 21 1272 

Wiese K J, 1998 Eur. Phys. ,J. B 1 269; ibid. 272 

Zimm B H, 1956 J. Chem. Phys. 24 269 

Klages R, Radons G and Sokolov I M (Eds.), 2008 

Anomalous Transport (Weinheim, Wiley- VCH) 

Scher H and MontroU E W, 1975 Phys. Rev. B 12 2455; 

Klafter J, Blumen A and Shlesinger M F, 1987 Phys. Rev. 

A 35 3081; Metzler R and Klafter J 2000 Phys. Rep. 339 

1 
[8] Mori H, 1965 Prog. Theor. Phys. 33 423; Kubo R, 1966 

Rep. Prog. Theor. Phys. 39 255; Zwanzig R, 1961 Lec- 
tures in Theoretical Physics (Boulder), Vol. HI, pp. 135 

(New York, Wiley) 
[9] Wang K G, 1992 Phys. Rev. A 45 833; Khan S and 

Reynolds A M, 2005 Physica A 350 183 

Lizana L et ai, 2010 Phys. Rev. E 81 051118 
[11] Metzler R and Klafter J, 2003 Biophys. J. 85 2776 
[12] Dubbeldam J L A et aZ., 2007 Phys. Rev. E 76 

010801 (R); Lua R C and Grosberg A Y, 2005 Phys. Rev. 

E 72 61918 

Chatelain C, Kantor Y and Kardar M, 2008 Phys. Rev. 

E 78 021129 

Zoia A, Rosso A and Majumdar S N, 2009 Phys. Rev. 

Lett. 102 120602 

Vocks H, Panja D and Barkema G T, 2009 J. Phys.: 

Condens. Matter 21 375105 
[16] Panja D, 2010 J. Stat. Mech. L02001 
[17] Panja D, 2010 J. Stat. Mech. P06011 
[18] Mandelbrot B B and van Ness J W, 1968 SIAM Rev. 10 

422; Lutz E, 2001 Phys. Rev. E 65 051106 
]19] Panja D and Barkema G T, 2009 J. Chem. Phys. 131 

154903 
120] Panja D, Barkema G T and Ball R C, 2007 



[10] 



[13] 

[14] 
[15] 



J. Phys.: Condens. Matter 19 432202; ibid. 

,arXiv:cond- mat/0610671, y2 
[21] Uhlenbeck G E and Ornstein L S, 1930 Phys. Rev. 36 

923 
[22] Kampen van N G, 2003 Stochastic processes in Physics 

and Chemistry (North-Holland, Amsterdam). 
[23] Rubinstein M, 1997 Phys. Rev. Lett. 59 1946 
[24] Heukelum van A and Barkema G T, 2003 J. Chem. Phys. 

119 8197; Heukelum van A et al., 2003 Macromolecules 

36 6662 
[25] Klein Wolterink J, Barkema G T and Cohen Stuart M 

A, 2005 Macromolecules 38 2009 



[26] Panja D. and Barkema G T, 2009 J. Chem. Phys. 131 

154903 
[27] Klein Wolterink J, Barkema G T, and Panja D, 2006 

Phys. Rev. Lett. 96 208301; Panja D, Barkema G T, and 

Ball R C, 2008 J. Phys.: Condens. Matter 20 075101; 

Panja D and Barkema G T, 2008 Biophys. J. 94 1630; 

Vocks H, Panja D, Barkema G T, and Ball R C, 2008 J. 

Phys.: Condens. Matter 20 095224 
[28] Panja D, Barkema G T, and Kolomeisky A B, 2009 J. 

Phys.: Condens. Matter 21 242101 



